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Abstract 

The conventional magnetic induction equation that governs hydromagnetic dynamo 
action is transformed into an equivalent integral equation system. An advantage of 
this approach is that the computational domain is restricted to the region occupied 
by the electrically conducting fluid and to its boundary. This integral equation 
approach is first employed to simulate kinematic dynamos excited by Beltrami-like 
flows in a finite cylinder. The impact of externally added layers around the cylinder 
on the onset of dynamo actions is investigated. Then it is applied to simulate dynamo 
experiments within cylindrical geometry including the "von Karman sodium" (VKS) 
experiment and the Riga dynamo experiment. A modified version of this approach 
is utilized to investigate magnetic induction effects under the influence of externally 
applied magnetic fields which is also important to measure the proximity of a given 
dynamo facility to the self-excitation threshold. 
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1 INTRODUCTION 



Dynamo action in moving electrically conducting fluids explains the existence 
of cosmic magnetic fields, including the fields of planets, stars, and galaxies pQ. 
As long as the magnetic field is weak and its influence on the velocity field is 
negligible we speak about the kinematic dynamo regime. When the magnetic 
field has gained higher amplitudes the velocity field will be modified, and the 
dynamo enters its saturation regime. 

The usual way to simulate dynamos numerically is based on the induction 
equation for the magnetic field B, 



where u is the given velocity field, fi the permeability of the fluid, and a its 
electrical conductivity. The behaviour of the magnetic field B in Eq. (CQ) is 
controlled by the ratio of field production and field dissipation, expressed by 
the magnetic Reynolds number R m = ficrLU, where L and U are typical length 
and velocity scales of the flow, respectively. When the magnetic Reynolds 
number reaches a critical value, henceforth denoted by R c m , the dynamo starts 
to operate. 

Equation ([T]) follows directly from pre-Maxwell's equations and Ohm's law in 
moving conductors. In order to make this equation solvable, boundary con- 
ditions of the magnetic field must be prescribed. In the case of vanishing 
excitations of the magnetic field from outside the considered finite region, the 
boundary condition of the magnetic field is given as follows: 



Kinematic dynamos are usually simulated in the framework of the differential 
equation approach by solving the induction equation For spherical dy- 
namos, as they occur in planets and stars, the problem of implementing the 
non-local boundary conditions for the magnetic field is easily solved by using 
decoupled boundary conditions for each degree of the spherical harmonics. For 
other than spherically shaped dynamos, in particular for galactic dynamos and 
some of the recent laboratory dynamos working in cylindrical geometry [2J, 
the handling of the non-local boundary conditions is a notorious problem. 

The simplest way to circumvent this problem is to replace the non-local bound- 
ary conditions by simplified local ones (so-called vertical field condition). This 
is often used in the simulation of galactic dynamos [3], and has also been 
tested in an approximate simulation of the Riga dynamo experiment [I]. 




V-B = 0, 





(2) 
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For the simulation of the cylindrical Karlsruhe dynamo experiment, the actual 
electrically conducting region was embedded into a sphere, and the region 
between the sphere and the surface of the dynamo was virtually filled by a 
medium of lower electrical conductivity [5115] . 



Of course, both methods are connected with losses of accuracy. In order to 
fully implement the nonlocal boundary condition, Maxwell's equations must 
be fulfilled in the exterior, too. This can be implemented in different ways. 
For the finite difference simulation of the Riga dynamo, the Laplace equation 
was solved (for each time-step) in the exteriour of the dynamo domain and 
the magnetic field solutions in the interiour and in the exterior were matched 
using interface conditions [7j. A similar method, although based on the finite 
element method, was presented by Guermond et al. [8119"] . Another, and quite 
elegant, technique to circumvent the solution in the exteriour was presented by 
Iskakov et al. [TOfTT] where a combination of a finite volume and a boundary 
element method was used to circumvent the discretization of the outer domain. 



An alternative to the differential equation approach (DEA) based on the so- 
lution of the induction equation is the integral equation approach (IEA) for 
kinematic dynamos which basically relies on the self-consistent treatment of 
Biot-Savart's law. For the case of a steady dynamo acting in infinite domains 
of homogeneous conductivity, the integral equation approach had already been 
employed by a few authors [T2|13|14|T5] . For the case of finite domains, the 
simple Biot-Savart equation has to be supplemented by a boundary integral 
equation for the electric potential [T6|rT7] . If the magnetic field becomes time- 
dependent, yet another equation for the magnetic vector potential has to be 
added [IE]. 



In the present work, the integral equation approach is applied to various dy- 
namo problems in cylindrical geometry. Two variants of the approach are 
presented: in the first one, it is implemented as an eigenvalue solver to solve 
genuine dynamo problems. In the second one, it is used to treat induction 
effects in the case of externally applied magnetic fields. Actually, the first 
variant was already at the root of the paper [19] where a surprising negative 
impact of sodium layers behind the propellers in the "von Karman sodium" 
(VKS) experiment was identified. It was not least this finding that prompted 
the VKS team to modify the experiment which made it ultimately succesful 
[20f2T] . After the derivation of the equation system in cylindrical geometry, 
we switch over to the treatment of specific problems, including the free decay 
case, the mentioned "von Karman sodium" (VKS) experiment [22]|23] . and the 
Riga dynamo experiment [241251126] . 



3 



2 MATHEMATICAL FORMULATION 



Assume the electrically conducting fluid be confined in a finite region V with 
boundary S, the exterior of this region filled by insulating material or vacuum. 
Then, dynamo and induction processes can be described [16] by the following 
integral equation system: 



v 



|kt f (u(r') x (Bo(r') + b(r'))) x (r - r>) ^ 



aaX r A(r') x (r — r') rTr . aa f , . .. . r — s' 



47T 



V s 



1 /■(u(r')x(B„(r')+b(r')) )_(s-rO (jr 



2 An J Is - r'| 3 



_ A / ^Hs-'V - JL / 0(sOn( S ') • ^> (4) 
47T 7 |s-r'| 3 An J |s - s'| 3 

A(r) = ±/W) + b(r-))x( r - r -)^, 
47r 7 |r - r'| 3 



where Bo is the externally applied magnetic field (which might be zero), b 
the induced magnetic field, u the velocity field, fi the permeability of the 
fluid (which is in most relevant cases the permeability of the vacuum), a the 
electrical conductivity, A the vector potential, and the electric potential, 
n denotes the outward directed unit vector at the boundary S. For a steady 
velocity field, the time dependence of all electromagnetic fields can be assumed 
to be ~ expAt. We have to distinguish three different cases: For non-zero Bo, 
and below the self-excitation threshold, the imaginary part of A is simply the 
angular frequency of the applied and also of the induced magnetic field. For 
B = the equation system (3-5) represents an eigenvalue equation for the 
unknown time constant A whose real part is the growth rate, and its imaginary 
part the angular frequency of the fields. For Bo = and A = 0, we need only 
the equations (3) and (4) which then represent an eigenvalue problem for the 
critical value of the velocity u at which the (non-oscillatory) dynamo starts 
to work. 
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2.1 General numerical scheme 



Although, in this paper, we will focus mainly on cylindrical systems it might 
be instructive to delineate the general numerical scheme for the solution of 
Eqs. (3-5). 

Assuming a specific discretization of all fields in Eqs. ([MS]), we obtain 



k = fia[P ik (B ok + b k ) - XRijAj - Qufa], (6) 
G m i4>i = S mk (B ok + bk) — XT m jAj, (7) 
A j = W jk (B ok + b k ), (8) 

where Einstein's summation convention is assumed. We have used the notion 
G m i = 0.5 o~ m i + U m [. B ok and b k denote the degrees of freedom of the externally 
added magnetic field and the induced magnetic field, Aj the degrees of freedom 
of the vector potential in the volume V, <pi the degrees of freedom of the electric 
potential at the boundary surface. Note that only the matrices Pi k and S mk 
depend on the velocity (the sources of the dynamo action), while Rij, Qu, T m j, 
G m \ and Wj k depend only on the geometry of the dynamo domain and the 
discretization details. 

Substituting Eqs.© and (jSJ) into Eq.flHJ) and eliminating Aj and <pi gives one 
single matrix equation for the induced magnetic field components bf 



h = fxa[P ik (B ok + b k ) — XRijWj k (B ok + b k ) — QuG lm S mk (B ok + b k ) 

+XQ u G br ]T m3 W jk (B ok + b k )} . (9) 

This equation can be further rewritten in the following form: 



[5 ik - [ioE ik - [iaXF ik \b k = [fiaE ik + fxaXF ik ]B ok , (10) 

where E ik = P ik - QuG^S mk and F ik = -R i:j W jk + QuG^T mj W jk . To com- 
pute induction effects, the induced magnetic field is obtained by solving the 
algebraic equation system (TlOl) . For the kinematic dynamo, Eq. lflO]) reduces 
to the following generalized eigenvalue problem 



[5 ik + [ioE ik \b k = X*F ik b k , (11) 
where A* is a new time constant rescaled according to A* = /icrA. 
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2.2 Cylindrical Geometry 



Since a number of dynamo experiments are carried out in cylindrical vessels, 
it is worth to specify the integral equation approach to this geometry. As long 
as the dynamo source (i.e. the velocity field or a corresponding mean-field 
quantity) is axisymmetric, the different azimuthal modes of the electromag- 
netic fields can be decoupled. This leads to a tremendous reduction of the 
numerical effort. The price we have to pay for this is the necessity to carefully 
deriving the dimensionally reduced version of the integral equation system. 

The electrically conducting fluid is assumed to be confined in a cylinder 
with radius R and height 2H. Introducing the cylindrical coordinate system 
(p, <f, z), we have 

r=[p cos <p, p sin <p, z} T , b = [b p , b v , b z ] T , u = [u p , u v , u z } T . (12) 

The magnetic field b, the electric potential <f), and the vector potential A are 
expanded into azimuthal modes: 



E 



m=—oo 



exp(im<£>). (13) 



When the velocity field is axisymmetric (i.e. it has only a component with 
m = 0), one can see that [b m , <f) m , A m ] T (m = 0,±1,±2, ••■) are decoupled 
with respect to m and they only depend on the variables (p,z). Henceforth, 
we always re-denote [b m , <fi m , A m ] T as [b, (f>, A] T for abbreviation. Then, after 
integrating over <p, Eq.Q acquires the form 



b p = ^[l \[{{z-zi)E™u z -ip'E™u v ){B 0p + b p ) + {-i{z-z')u z E, 



rn 



H R 
-H 

+ip'E™u p )(B 0v + b v ) + (i(z - z')E^ -(z- z')u p E™)(B 0z + b z ))p'dp'dz' 

R H R 

- J <j)p' 2 E™\ z , =H dp' - J <j>R{z-z')E^\ pl=R dz' + J <j> P a E™\ z , = _ H d P ' 

-H 

H R 

-A f ({{z-z')E™A p + {z-z')E™A v + p'E™A z )p'dp'dz% (14) 



-H 
H R 

K = ^iJ J [{-{PK ~ P'KW + i(* - z')u z ET)(B 0p + b p ) 

-H 

+ (( P E 1 m - p'E?)u p + (z- z')u z E™){B 0tp + b v ) + {-(z - z')u v E r c 
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R 



-i(z - z')u p E™)(B 0z + b z )]p'dp'dz' - J <f)(pp'E™\ z , =H - p l2 E™\ z , =H )dp' 

o 

H R 

+ J <t>R{z-z')E™\ p , =R dz'- f<i>(-pEr\ z , = -. H + p , E?\ g , = - H )f/df/ 

-H 
H R 

-X f [ {{z - z')E™A 9 -{z- z')E r c n A p + {pE? - p 'E™)A z )p 'dp 'dz'}, (15) 



-H 
H R 



bz= 4^ [ J I [{p>K ~ pE ^ Uz ( B °» + ^ + wKuziB^ + b v ) 

-H 

H 

+ ((- p 'E? + pE™)u p -ipE™u lp )(B Qz + b z ))p'd P 'dz' + J 4>RpET\ P >=Rdz' 

-H 

H R 

J rpm „jpm\ a \ J J J J ~'l 



-A j ]{-pE™A p +{p'E™-pE™)A 9 )p'dp'dz'l (16) 

—H 

where the following azimuthal integrals appear: 



27T 

P m/ / a / cosm^ , 

^ (p,p ,z,z )— / 3-a</9 , 

J (p 2 + p' 2 — 2pp' cos ip' + (z — z') 2 )2 

E™(p,p',z,z')= / " w dip', 

J (p 2 + p' 2 — 2pp' cos y?' + (z — z') 2 ) 2 

i£ s (p,p ,z,z)= / jcty? . 

J ( n 2 -I- n' 2 — 9nn' ms in 1 A- (7. — 7. 1 ) 2 ) 2 



Accordingly, from Eq.(jH), we obtain the expressions for the electric potentials 
at the three different surface parts of the cylinder: 



1 h R 

2 0(Sl) = 4^ [ / / i-p'P E ?\z= H u z - p'{H - z') Utp E?\ z=H )(B 0p + b p ) 

-H 



+((-p'pE™\ z=H + p' 2 E™\ z=H )u z + p\H - z>)u p E™\ z=H )(B 0v + b v ) 



— / (j)R(pE™\ p > =R>z=H — RE™\pi =RjZ=H )dz' 

-H 

R H R 

+2.0H J <pE™\ z=H , z , = _ H p'dp' - A J J p'(pE™\ z=H - p'E r {% =H )A p 



-H 
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-p'pE™\ z=H A v + p'(H - z')E™\ z = H A z dp'dz'\, (17) 



1., , 1 H R 



2^) = ^[j I (-p'RE^\ p=R u z - p'(z - z')u v E?\ p=R )(B 0p + b p) 



-H 

„/2 771m I 



+(-p'RE™\ p=R u z + p' 2 E™\ p=R u z + p'(z - z')u p E™\ p=R )(B 0v + b v ) 
+ (p'u v RE™\ p=R + p'Ru p E™\ p=R - p' 2 u^E™\ p=R ){B Qz + b z )dp'dz' 

R H 

2jJ 



~ j <p(z-H)ET\ P =R, z '=Hp'dp'~ J <j)(E™\ p=p , =R -E™\ p=p , =R )R 2 dz 

—H 
R H R 

+ J<j ) (z + H)E™\ p=R , zl= _ H p'dp'-\ J J(p'RE™\ p=R 



-HO 
J2 rpm I 



-p a E?\ p=R )A p - p'RE™\ p=R A v + p\z - z')E™\ p=R A z dp'dz\ (1* 



1 " />> 



2^ s 3) = ^[ / / (-p'p^ m |,=-^ + p'(# + z')ET\ z =-hu v )(B 0p + b p 



-H 

12 rpm I 



+ (-p'p^c m |,=-^ + P 2 ET\ z =-hu z - p'(H + z')u p E™\ z= _ H )(B 0ip + b v ) 
+ (p'pu v E™\ z= _ H - p' 2 u v Er\ z =-H + p' P u p ET\ z =-h)(B 0z + b z )dp'dz' 

R H 

+2.0H J <pE™\ z= _ HyZl=H p'dp' - J <j)R{pE7?\ p , =RjZ= ^ H 



-H 



-RE™\ p , =R , z= _ H )dz' - A J {pp'E™\ z= _ H - p' 2 E?\ z= _ H )A p 

v 

-p'pE™\ z= _ H A v + p\-H - z')A z ET\ z= - H dp'dz'}. (19) 



Here, s\ is the surface z = H, s 2 the surface p = R, s% the surface z = —H. 
Equation (jSJ) for the vector potential gets the form 



1 " 



A p = ^U I P'( z - z ')K{Bo P + b p ) + p\z - z')E™(B 0ip + \ 



-H 

R 



+p' 2 E^(B 0z + b z )dp'dz' + J -p'DT\ z , =H (B 0p + b p ) - p'D™\ z , =H {B Q(p + b v )dpl 

o 

H R 

+ J R(B 0z + b z )D™\ p ,= R dz' + J p'DT\ z '=- H (B 0p + b p ) 

—H 

+p'(B 0v + b v )D™\ z , = _ H dp'] (20) 

K = ^\\ J -A* - + b p ) + p\z - /)^r(^ + K) 

-H 
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+p'(pE™ - p'E™)(B 0z + b z )dp'dz' + J p'D™\ z , =H (B 0p + b p ) 

o 

H 

-p'D™\ z , =H (B^ + b^)dp' - J R(Bv z + b z )D™\ p>=R dz' 

-H 

R 

+ J p'D™\ z , = _ H (B 0v + b v ) - p'D™\ z , = _ H (B 0p + b p )dp'\ 
o 

H R 

Az = ^ [ I I -P'P E ™( B °P + h p) + P'(P' E ™ ~ P E c)(B 0lp + b v )dp'dz' 

-H 
H 

+ J RD™\ pl=R {B^ + b^)dz% 

-H 

where the following abbreviations of azimuthal integrals were used: 



2tt 



D: i (p,p',z,z') 
D™{p,p',z,z') 
D?(p,p',z,z>) 



sin <p' sin mip' 



J (p 2 — 2pp' cos <p' + p' 2 + (z — z') 



-dip' 



2n 



cos ip cos rmp 



J (p 2 - 2pp' cos <f' + p' 2 + ( 



z — z') 2 ) 2 



dtp' 



2n 



cos my? 



•> (p 2 — 2pp' cos ip' + p' 2 + (z — z') 



-dip 1 . 



In our numerical scheme, we typically use equidistant grid points pi = i x Ar 
and Zj = j x A;z to discretize the intervals [0, i?] and [-H, H], respectively (in 
some applications non-equidistant grid points are also used). The extended 
trapezoidal rule is applied to approximate all the integrals in Eqs. (I14ti22p . 
Then we obtain the following matrix equations 



1 

2 



f A 

A 

\A Z 







f B 0p + b p \ 


(<t>sl\ 




(A\] 


pa 


P 


Bq v + b v 


- Q U»2 


- AR 


A 






\ B 0z + b z ) 






\A)_ 



Bq p + b p 
Bop + b<p 
B 0z + K 
Bq p + b p 
W | 1'> {W + b. 
Bqz + b z 



A 

AT I 

A, 



U 



/ <f>8l 

i>s2 



(23) 
(24) 
(25) 
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where the matrix elements of P, Q, R, S, T, U, and W can be read off from 
Eqs.(H3H22P. Combining Eqs.fl23H25J, we obtain 



(I - fwE - na\F)b = //tr(E + AF)B, 



(26) 



where 




(27) 



(28) 



After solving the algebraic equation system (1261) . the induced magnetic field 
b can be obtained for the magnetic induction process. 

For the kinematic dynamo problem, the following generalized eigenvalue prob- 
lem has to be solved 



for the given velocity field u, where A* = fiaX. Note that a quite similar 
numerical scheme can be established in spherical geometry for the case of 
axisymmetric dynamo sources. 



3 NUMERICAL IMPLEMENTATION AND RESULTS 

In this section, the integral equation approach will be applied to various cylin- 
drical dynamo and induction problems of experimental relevance. We start 
with the problem of the free decay of a magnetic field in a cylinder. Then, a 
class of Beltrami-like flows will be considered. In all the problems we use the 
QZ algorithm [27] which is a modification of the QR algorithm for the case of 
generalized non-hermitian eigenvalue problems. 

The integral equation approach is further employed to investigate the induc- 
tion effect of the VKS experiment. The algebraic equation system is solved by 
the LU decomposition. The obtained induced magnetic field will be compared 
with the data measured in experiment. 

At the end we deal with the Riga dynamo experiment with its large ratio of 
height to radius. Due to the large resulting matrices we shift here from direct 
matrix inversion methods to the generalized inverse iteration method |28j . 



(I - fwE) ■ b = A*F ■ b 



(29) 
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Fig. 1. Freely decaying magnetic field in a finite cylinder with R = 2H = 1. 

3. 1 Free field decay in a finite cylinder 

The simplest problem to start with is the free decay of a magnetic field in a 
finite length cylinder. This example was already treated by Iskakov et al. |10j . 
In Fig. [1] we show the magnetic field lines of the slowest decaying eigenfield, 
which has the same dipolar structure as in Fig. 8 in [ID]. 

3.2 Beltr ami-like flows 

In this subsection, we consider a class of flows in finite cylinders which we call 
" Beltrami-like" flows. Actually, the velocity field v(r) of Beltrami flows are 
characterized by the property V x v = f3v. In dynamo theory, there is par- 
ticular interest in such flows since they are also helicity maximizing. Helicity 
maximizing flows, in turn, are well known to possess quite small critical R m , a 
fact that was utilized, e.g., in the optimization of the Riga dynamo experiment 
[7]. The actual flow structures that will be treated in this work were proposed 
by J. Leorat [29] , and a certain sub-class of them (with ideally conducting 
boundary conditions, however) was considered by Wang et al. [3D] . 

We use the notation s^t n to characterize flows with m poloidal vortices and n 
toroidal vortices. The sign ± indicates that the poloidal flow in the equatorial 
plane is directed inward (+) or outward (— ), respectively. An impression of 
the topology of the flow structure can be obtained from Fig. [2] where we have 
also indicated possible propeller or rotating disk configurations to produce 
such flows. 

The analytical expression of the flows to be considered in this paper are as 
follows: 
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Fig. 2. Illustration of the considered flow topologies s^t n , and of typical propeller 
or rotating disk configurations to produce them. 



v r (r, z) = c\ J±(a r) cos(m7r(,2 + H)/2H), 

v v (r, z)=t Ji(a r) cos{nTt(z + H)/2H), 

v z (r, z) = —c\ c 2 a /ir J (a r) sm(mn(z + H)/2H), 



(30) 
(31) 
(32) 



where a = 3.8317 is the first root of the Bessel function J\. In the following, we 
will restrict ourselves to m,n = 1,2. c\ — 1 for all s+i n flows, c\ = —1 for all 
s~t n flows, C2 = H/2 for s\t\ flow and for other flows Ci — 1. Again, 2if is the 
height of the cylinder, 2; G [— H, H] . In the following discussion the height of 
the cylinder is set to 2 and the radius is fixed to 1. The parameter r indicates 
the ratio of toroidal to poloidal flow. We will consider values of r close to 2 
which turned out to be advantageous for dynamo action. Actually, this value 
r = 2 is not the value which would correspond to an excact Beltrami flow. 
This is the reason why we have called the considered flows " Beltrami-like" . 

In what follows we will use a definition of the magnetic Reynolds number R m 
which is based on the maximum of the axial velocity. In order to display the 
results for the s^t n in one common figure we will use the ± as a sign of R m 
according to i? m = ±/iO"i?|u^ a3; |. 

For the case without external layer (w = 0) we found that only the 
dynamo is steady, all the others are all oscillatory. However, if an external 
layer around the finite cylinder is added, even the thickness of the layer is 
quite small, for example, equal to 0.05, the s2~t2 dynamo becomes steady. 

The magnetic field structures for sjii, s^ti, s^h, sJT^ and flows are 
shown, at an azimuthal section at cp — 0, in Figs. El El El El and El, respectively. 
In all these cases, an externally added layer with thickness equal to 0.5 has 
been considered. The variations of growth rates of the magnetic fields with 
respect to the magnetic Reynolds number for all the flows are depicted in 
Figs. HI [7| and [Hfl From these figures, one can see that the externally added 
layer has a very strong impact on the onset of dynamo actions. For example, 
the critical magnetic Reynolds number for the flow s^ti is approximately 
equal to 143 in the case without external layer. When an external layer with 
thickness 0.2 is considered, the critical magnetic Reynolds number reduces 
to 61. If the thickness of the external layer is increased to 0.5, the critical 
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15 




Fig. 4. Growth rates for s±ti flow and the influence of the externally added layers 
with thickness w. 

magnetic Reynolds number further declines to 40. Finally, one can also note 
that there is a tendency that when the thickness of the external layer becomes 
larger, the curves of the growth rates become more symmetric with respect to 
the ordinate axis at R m = 0. 

In Fig. [IT] we show, for the flow with different widths w of the external 
layer, the dependence of the critical R m on the parameter r which measures 
the ratio of toroidal to poloidal motion. For w = we show, in addition to the 
results of the integral equation approach, also the results of a finite difference 
code based on the differential equation approach as it was described in [7] and 
also used in [19]. In general, we observe a good correspondence of the results 
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3.3 Induction effects in the VKS Experiment 

Since the flows are characterized by a comparatively low critical R m , 
much focus was laid on their realization in experimental dynamos. Both the 
spherical Madison dynamo experiment (MDX) p?Tll52] and the cylindrical "von 
Karman sodium" experiment (VKS) [22.23J are realizations of the s^tz flow. 

Although a recent version (using impellers with a high magnetic permeability 
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Fig. 8. Magnetic field structure for s 2 £2 flow 



fi re i ~ 200) of the VKS experiment has shown dynamo action [20] and even 
a kind of polarity reversals [21], the under-performance of the original VKS 
experiment compared to numerical predictions is still a matter of interest. In 
[T9] we have tried to explain this by the detrimental effect of sodium layers 
behind the propellers on the dynamo action. The sheer existence of these 
layers leads already to a significant increase of the critical R m which becomes 
dramatic if a realistic rotation therein is taken into account. 



Interestingly, in the original VKS experiment, the measured induced magnetic 
fields, for large R m , are significantly weaker than the numerically predicted 
ones. Using our method we will try to figure out if this effect can also be 
attributed to the existence of lid layers and the flow therein. 



15 



(a) 



(b) 



(c) 



Fig. 9. Magnetic field structure for s^t? flow. 
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Fig. 10. Growth rates for s 2 £2 flow an d influence of the externally added layers 



A rather realistic flow field, resulting from the so-called TM73 propeller [34 
(which was identified as a sort of optimal flow field), is considered in our 
calculation. Some interpolations were necessary to project this flow field onto 
the grids used in our code. More details on this can be found in [19|. 

In the original VKS experiment a static external magnetic field was applied 
in the transverse direction. The induced magnetic field was measured in the 
direction perpendicular to the externally added magnetic field at the point r = 
0.5 in the equatorial plane. In the following the induced magnetic fields near 
the r = 0.5 points obtained by our integral equation approach are compared 
with the measured ones. The influence of a rotating flow in the lid layer on 
the induced field is investigated. Three kinds of velocity field in the lid layer 
are considered. The first one is a static lid layer. The second one is that only 
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Fig. 11. Critical R m for the s^~t2 How in dependence on the toroidal/poloidal ratio 
r for different thicknesses w of the surrounding layer. For w = 0, we compare the 
results with the results of a differential equation solver. 
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Fig. 12. Ratio of induced axial (left) and azimuthal (right) magnetic field to the 
applied magnetic field. Experimental results at radius r = 0.5 (taken from [35]) and 
numerical results under the assumptions of different flows in the lid layer. Static 
lid layer, constant velocity in lid layer, constant velocity in lid layers multiplied by 
factor 1.5. 

a rotation of the lid layer is assumed, but v v remains constant in the axial 
direction, its dependence on the radial variable is the same as on the interface 
between the lid layer and inner part of the cylinder. The third one has only 
one difference from the second case in that the magnitude of v v is increased 
by a factor 1.5 which comes closer to the velocity of the impeller. For these 
three cases, the numerical axial induced fields around the point r = 0.5 and 
the experimental result at r = 0.5 on the equatorial plane are shown in the left 
part of Fig. [12] for different magnetic Reynolds numbers. The right panel of 
Fig. [12] displays the azimuthal induced magnetic fields. From these figures one 
can see that the third case shows a best agreement with the experimental one. 
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(a) (b) 

Fig. 13. Simulated magnetic field structure of the eigenmode of the Riga dynamo 
experiment, (a) Isosurface plot of the magnetic field energy. The isosurface corre- 
sponds to 25 per cent of the maximum magnetic energy, (b) Magnetic field lines. 

A good agreement of the axial magnetic field with the experimental result has 
been achieved for the second and third cases. But for the azimuthal magnetic 
field, when the magnetic Reynolds number is larger than 30, there is still a 
gap between the numerical results and the experimental ones. 

Nevertheless it is quite likely that it is indeed the existence of lid layers and 
some azimuthal flow therein which is responsible for the unexpected under- 
performance of the original VKS dynamo experiment. 

3.4 Riga experiment 

In this subsection, the integral equation approach is used to re-simulate the 
kinematic regime of the Riga dynamo experiment. This experiment has been 
optimized and analyzed extensively within the differential equation approach 
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Fig. 14. Comparison of the IEA and DEA results for the Riga dynamo experiment, 
together with experimental results, (a) Growth rate, (b) Frequency. 

(DEA) by means of a finite difference solver [7]. The values of the velocity field 
on the grids used in our code are obtained by interpolating the experimental 
velocity field measured in a water-dummy experiment. The influence of less 
conducting stainless steel walls has not been taken into account. 

The computations have been carried out on a 100x20 grid in z- and r-direction. 
The structure of the magnetic eigenfield is illustrated in Fig. [T5J Figure [T3"a 
shows the isosurface of the magnetic field energy (this time at 25 percent of the 
maximum value). In Fig. [13b the magnetic eigenfield lines are depicted. Basi- 
cally the structure is the same as that resulting from the differential equation 
approach [26] with a 401x64 grid in z- and r-direction. 

The dependence of the growth rate and frequency of the eigenmode of the Riga 
dynamo experiment on the rotation rate is shown in Fig. [T4"r and Figfl4"b. 
respectively. The comparison with the DEA results shows that the slopes 
of the curves are in good agreement. However, we see that the limited grid 
resolution in the IEA leads to significant shifts in the order of 5 percent towards 
lower rotation rates for the growth rate and of 10 per cent towards higher 
rotation rate for the frequency. Hence, it could be said that the Riga dynamo 
experiment marks a margin of reasonable applicability of the IEA with its 
need to invert large matrices which are fully occupied. 



4 CONCLUDING REMARKS 



In the present paper, the integral equation approach to kinematic dynamos has 
been applied to non-spherical geometries. The method was examined by its 
application to the free field decay. The comparison of the obtained results with 
other methods shows a good agreement. The integral equation approach was 
extended to investigate induction effects of the VKS experiment. The obtained 
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induced magnetic field shows a satisfactory agreement with the experimental 
result when the effect of the lid layers and a certain azimuthal flow therein 
are taken into account. Finally, it was applied to simulate the Riga dynamo 
experiment. 

It can be concluded that the integral equation approach is robust and reliable 
and can be used for practical purposes, although limits of its applicability are 
seen for the Riga dynamo experiment with its large ratio of length to radius. 
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